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A body immersed in a highly viscous fluid can locomote by drawing in and expelling fluid through 
pores at its surface. We consider this mechanism of jet propulsion without inertia in the case of 
spheroidal bodies, and derive both the swimming velocity and the hydrodynamic efflciency. Ele- 
mentary examples are presented, and exact axisymmetric solutions for spherical, prolate spheroidal, 
("^^ and oblate spheroidal body shapes are provided. In each case, entirely and partially porous (i.e. 

T-H jetting) surfaces are considered, and the optimal jetting flow profiles at the surface for maximizing 

^—^ the hydrodynamic efficiency are determined computationally. The maximal efficiency which may be 

achieved by a sphere using such jet propulsion is 12.5%, a significant improvement upon traditional 
>^ > flagella-based means of locomotion at zero Reynolds number. Unlike other swimming mechanisms 

^^ which rely on the presentation of a small cross section in the direction of motion, the efficiency of 

^-i a jetting body at low Reynolds number increases as the body becomes more oblate, and limits to 

^^ approximately 162% in the case of a flat plate swimming along its axis of symmetry. Our results 

are discussed in the light of slime extrusion mechanisms occurring in many cyanobacteria. 
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I. INTRODUCTION 

Locomotion at the micron scale, or in highly viscous fluids, is constrained by the dominance of viscous dissipation 
over inertial effects. Hence, the strategies utilized by microorganisms and engineered swimming devices must differ 
from more familiar locomotive mechanisms such as the flapping of wings [1-4] . Instead of imparting momentum into a 
fluid wake, which is not possible at zero Reynolds number, bodies generally exploit drag anisotropy in order to propel 
themselves through the fluid. In particular, the undulation of flagella and cilia are the most well-studied means of 
swimming at low Reynolds numbers in nature. As shown in classical explorations of flagellar locomotion [5, 6], the 
maximum theoretical hydrodynamic efficiency which may be achieved by flagellar propulsion is approximately 8%, 
and is generally closer to 2 % in real cells [5, 7, 8]. For locomotion using cyclic tangential surface distortions, the 
efficiency has a theoretical bound of 75% in the case of a spherical body, a dramatic improvement upon that granted 
by the use of flagella [9] . 

Flagellar locomotion is not the only means of locomotion at low Reynolds numbers and many authors have considered 
alternatives, either to help explain biological phenomena, or to suggest designs for synthetic locomotor systems on 
small scales. Examples of swimming bodies which deform in a manner which breaks the important time-reversal 
symmetry (a constraint known as the Scallop theorem and without which no locomotion is possible in the absence 
of inertia) were presented by Purcell [2]. These systems include the motion of a three-link swimmer [2, 10-12] and 
a treadmilling torus [2, 13-16]. Recently, Leshansky et al. [17] have considered the remarkably efficient locomotion 
of an elongated treadmilling body, which can propel itself at nearly the same velocity of the surface motion. Other 
simple means of propulsion proposed for locomotion at zero Reynolds number include a deformable two-dimensional 
loop [18], systems of two [19, 20], three [21, 22] or N spheres [23], flapping near deformable interfaces [24], a rehinging 
swimmer [25], and a jellyfish-inspired bilayer vesicle [26]. For a more complete list of references we refer the reader to 
Ref. [7]. 

Recently, it has been observed that certain types of bacteria, such as myxobacteria and cyanobacteria (blue-green 
algae) secrete mucilage through nozzle-like organelles while gliding along a substrate (see Fig. 1) [27-32]. This so- 
called "slime extrusion" has been theorized as a primary propulsion mechanism for adventurous motility in such 
organisms as M. xanthus. The slime is a polyelectrolyte gel, and has been modeled after snail slime. Wolgemuth et 
al. [31] have shown that the osmotic expansion of the slime from the nozzle generates a sufficient force to propel the 
organism. As the slime exits the nozzle, it adheres to the substratum, and further slime extrusion produces a thrust. 
Though the primary activation is located near the poles of the elongated body, the entire surface is covered with 
nozzles. Meanwhile, the motility mechanism employed by the motile marine cyanobacteria Synechococcus is still an 
open problem in biophysics [9, 33-35]. Synechococcus is known to swim absent the presence of a substrate, without 
changing shape, and without any observable external organelles, a puzzle to which we shall return. 

Inspired by the gliding locomotion of the cyanobacteria, we consider in this paper the following questions: Can the 
extrusion of a Newtonian fluid be used as a propulsive mechanism absent a solid surface? And if so, how efficient 
would such a swimmer be? Herein we study the swimming motion which may be achieved by the placement of surface 
nozzles upon a body surface which act to both draw in and expel fluid, as illustrated in Fig. 2. We refer to this 
locomotion as extrusion swimming, or non-inertial jet propulsion. The swimming mechanism is not unlike the more 
familiar jet propulsion at higher Reynolds numbers, but with distinctly different flow structure and resultant (non- 
local) fluid-body interactions. In classical high Reynolds number jet propulsion, a body propels itself by imparting 
momentum onto the fluid opposite the direction of motion (as in the swimming of jellyfish) [4, 36-42]. The low 
Reynolds number analogue studied here propels itself instead by taking advantage of the viscous stresses induced by 
the jet motion on its surface. 

The fluid being expelled from the body is assumed to be Newtonian and identical to the surrounding fluid, and 
the body is assumed to be well-separated from any surfaces. The swimming velocity and hydrodynamic efficiency of 
such a body arc first derived formally using the Lorentz reciprocal identity. We find that the maximal hydrodynamic 
efficiency which may be achieved by a sphere using such jets is 12.5%, a significant improvement upon most other 
means of swimming at zero Reynolds number. Moreover, unlike many other swimming mechanisms which rely on 
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FIG. 1. Reprinted from Wolgemuth et al. [31] with permission from Elsevier, (a) Negatively stained electron micrograph of 
an isolated M. xanthus cell envelope showing multiple ring-like structures located predominantly at the poles of the cell. The 
inset shows a higher magnification of the nozzle array in the region indicated by the long arrow, (b) A gallery of electron 
micrographs of negatively stained isolated nozzles. Each cylindrically symmetric nozzle has an outer diameter of 14 nm. (c) 
Electron micrograph of a gliding cell. At higher magnification, it can be seen that the slime trails are composed of several slime 
bands, which are secreted from the sites at the cell pole, where the nozzles are located (large arrow). 



the presentation of a small surface area in the direction of motion, we show that the efficiency increases as the body 
becomes more oblate, and limits to approximately 162% in the case of a flat plate moving along its axis of symmetry. 
The manuscript is organized as follows. In §11, the swimming velocity of an arbitrary body is shown to be dependent 
upon the fluid stress in a dual (but simpler) resistance problem, and we present elementary examples for certain 
swimming spheroids and a rotating "viscous pinwheel." The hydrodynamic efficiency of locomotion is also defined. 
In §111, we consider the behavior and efficiency of a spherical body, for the case of an entirely porous (i.e. jetting) 
surface, and then for a partially porous surface. We also determine the optimal jetting flow profile at the surface for 
maximizing the swimming efficiency. In §IV and §V we consider prolate and oblate spheroidal bodies, respectively. 
Entirely and partially porous surfaces are considered, and jetting profile optimization is performed. We conclude by 
applying the theoretical results to the organism Synechococcus, which swims by a yet-unknown mechanism, and show 
that while slime extrusion is important for gliding motility in related cyanobacteria, extrusion of a Newtonian fluid 
cannot not properly account for its observed swimming speed. 



II. GENERAL EXPRESSIONS FOR THE SWIMMING VELOCITY AND EFFICIENCY 
A. Description of fluid-jetting bodies and fluid-body interactions 



We begin our consideration by calculating the swimming velocity of a jetting body of arbitrary shape. A fluid- 
jetting body is illustrated in Fig. 2. The surface is composed of a solid part, ODq, and the m porous surfaces which 
allow for inflow and outflow, dDg for i — 1,2, ...,m. A no-slip condition is assumed to hold on the surface dD^), and 
the porous flow through the surfaces dD\ is assumed to be driven by internal mechanisms in the direction normal to 
the surface at each point, either into or out of the body. The entire surface is denoted by dD = dDo U™ j^ dDg. 

The equations for incompressible fluid motion at zero Reynolds number are the Stokes equations. 



V -cr = -Vp + /iAu = 0, 
V • u = 0, 



(1) 
(2) 



where cr — —pI+2^E is the Newtonian fluid stress, p is the dynamic pressure, u is the fiuid velocity, E — (Vu+Vu-^)/2 
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FIG. 2. (a) A fluid-jetting body. The surface is composed of a solid part, 9Do, and the m porous surfaces which aUow for 
inflow and outflow, dD\ for i = 1, 2, ..., tti. The entire surface is denoted by dD — dDo U"Li dDg. Here the body draws in and 
expels fluid to the left, and swims to the right with velocity U. We denote by n the unit normal vector pointing into the fluid. 
(b) The surface dD} corresponds to the i"* fluid jet, or porous surface. 



is the symmetric rate-of-strain tensor, and fi is the shear viscosity. The boundary conditions are decomposed into 
rigid body motion and a fluid extrusion (or jetting flow) component. 



u(x) = U + f2 X X + </«(x)n(x) (xeaD), 
u(x) — > 0, p -^ (x-^oo). 



(3) 
(4) 



The jetting-flow profile 0(x) has support only at the porous surfaces, x G U™ ;^9£'J, or 0(x) = for x e dDo- We 
denote by n(x) the unit normal vector at the point x pointing into the fluid, and U and Jl are the instantaneous rigid 
body velocity and rotation rate resulting from the fluid extrusion. We define the fluid-jetting fluxes (i.e. flow rates) 
Qi via 



/ (t){x)dS^qi, 
JdDi 

with dS the surface area element. The internal volume is conserved if the inward and outward fluxes balance, 

/ n.(</.(x)fi)d5 = Vg, -0. 

Finally, to close this "swimming problem" we require that the body be force and torque free, i.e. 

/ (TxidS = 0, I xx{cT -n) dS = 0. 

JdD JdD 



(5) 



(6) 



(7) 



B. General expression for the swimming velocity 



An identity attributed to Lorentz may be used to deduce the rigid body motion resulting from the porous jetting 
flow on the body surface. This approach was used by Brenner to compute the drag on a body in an arbitrary flow 
field [43] , and by Stone and Samuel to study the swimming of finite size bodies [9] . The Lorenz reciprocal theorem is 
stated as 



/ u ■ {cr ■ h) dS = u • (cr • n 

JdD JdD 



dS, 



(8) 



where (u, cr) are the velocity and stress fields corresponding to the swimming problem described above, and (u, a) 
are the velocity and stress fields corresponding to a different system, but one which shares the same instantaneous 



immersed boundary (see Refs. [44, 45]). Let (u, a) be the solution to the Stokes equations for rigid body motion of the 
immersed body, so that the surface fluid velocity isu = U + i7xx and the corresponding net forces and torques are 
F and L respectively (the "resistance problem"). Since the body is force and torque free in the swimming problem, 
the integral on the left hand side of Eq. (8) vanishes. The remaining terms in the reciprocal identity are 

0=/ (U + fixx + (f)(x) n) • (ct . fi) dS. (9) 

JdD 

Since 0(x) = for x e dDo, we obtain 

U-F + rj-L = - / 0(x) {h ■ a ■ h) dS = -V] / (/)(x) (n • 5r . n) dS. (10) 

JdD i^i JdDi 

In the event that a body is to translate or rotate with the greatest velocity for a given set of jet fluxes qi, Eq. (10) 
indicates that the jets should be placed precisely where the body experiences its largest normal component of traction 
in the corresponding resistance problem. Also of note, the consequences on the flow and swimming velocities from 
inward flowing jets are precisely opposite those of outward flowing jets. Unlike at higher Reynolds numbers where 
there is a distinct asymmetry in the flow fields set up by an outward flux and inward flux through a small opening 
(candles are blown out instead of sucked out!), at zero Reynolds number sources and sinks near a wall produce identical 
streamlines [44]. We proceed to consider a number of simple examples which utilize Eq. (10). 



C. Example 1: A translating sphere 

As a first example, consider the simple case of a translating sphere of radius a. The corresponding resistance 
problem is well known: F ~ 6TT^a\J, L = 0, and {a ■ h) ~ 3/i/(2a)U on the body surface [44]. Generally, then, the 
swimming velocity of a jetting sphere may be written as 

V = -^J^H^)ndS. (11) 

The largest swimming velocity for a balanced inward/outward fluid flux is therefore achieved by placing two jets of 
equal and opposite strengths at the spherical poles along the axis of locomotion (xi and X2), 

0(x) = rj((5(x,xi)-<5(x,X2)), (12) 

where (5(x, y) is the Dirac delta function with support at x = y, and q = qi = —52- The corresponding swimming 
velocity is 

While the flow profile 4>(x) stated above yields the largest swimming velocity, it does not present the most efficient 
means of moving the sphere through the fluid, as we shall explore in §111. 

More generally, if the pores are small, Sg -^ 1 (with S} the surface area of the i*'* porous surface dDg), then it is 
useful to expand the stress from the resistance problem linearly about the jet locations. For instance, near the i*'* 
porous jet we have 

n(x) = fi(x,)+0(^l), (14) 

(^•n)(x) = (5-.n)(x,) + 0(5^), (15) 



and Eq. (10) may be written as 

m 

V-F + n-L = -y^qi (n-CT-n)(x,) + 0(max5'j). (16) 

The swimming velocity of a sphere with m smah fluid jets acting upon its surface is, from Eq. (16), 

i=l 

A jetting sphere may therefore move in any direction spanned by the normal vectors at the jet locations by tuning 
the fluxes qi, and the resulting motion will be a simple translation in that direction. Regardless of their distribution 
and strengths, jets acting normal to the body surface cannot, however, be used to generate rotations of a sphere. For 
the case of a rotating sphere of radius a, the resistance problem has F = 0, L = Siriia^Ct, but also (n • o- • n) = 0, so 
that Eq. (10) gives H = regardless of the jetting flow profile 0(x). In order to rotate using porous extrusion, the 
body must not be axisymmetric about the axis of rotation. 



D. Example 2: A spheroid translating along its axis of symmetry 

The swimming velocity of a jetting spheroid translating along its axis of symmetry can also be determined with 
ease. Once again the corresponding resistance problem has a long history [44]. Consider a prolate spheroid with major 
and minor axis lengths 2a and 26, with its major axis aligned with the z axis. Setting U = Ui, we have F = 67r/i_RU, 
L = 0, and 



2^C/C 



Vc(r2-C2)[(r2 + l)coth-i(ro)-ro]; 

on the body surface, where C = cos(0) [6 is the polar angle), c — vc? — 1?, tq — a/c, and 

8c/3 



(18) 



R 



(t2 + 1) log 



To + 1 



(19) 



-2to 



.^0 - 1. 
(see [44]). Equation (10) then gives the swimming velocity; after some algebra we find that 



1 -^2_A2 



1 ^^..^...^.. ^..^_../l_-C^ 



U = '2 I ^(^)'^(^) '^^' ^(^) = ^ V 7^ • (20) 



The function G(C) is monotonically increasing in ^ for all values of tq with G(G) = and G(l) = 1. Confirming 
intuition, jets placed nearer to the poles contribute more significantly to the swimming velocity. Placing jets only at 
the poles, which expel and draw in fluid with fluxes ±q, the swimming velocity (with U = Uz) is 

For b = a we recover the spherical swimming velocity, and as the body becomes more slender the swimming velocity 
increases without bound (for fixed q). This is the largest velocity which may be achieved by a jetting prolate ellipsoid 
in the direction of its major axis for a given inward/outward flux, and as we will show in §IV is again not the most 
efficient. 

The swimming velocity of an oblate ellipsoid translating along the symmetric axis can be obtained by applying to 
the above the transformation (c, tq) —>■ (i c, —i Aq), and reversing the definitions of a and b so that Aq = b/y/a'^ — 6^ > 



(see Ref. [44]). In this case the body swims with velocity 



U = 



Ci 



'A§ + C^ 



A^ 



1 



0(C) dC, 



(22) 



which is again maximized by placing jets at the poles passing fluid with fluxes ±g, giving U — —q/{2iTa'^). For fixed g, 
the swimming velocity decreases without bound as the presented surface area and hence fluid drag (due to the large 
no-slip surface area) increases. 



E. Example 3: A prolate spheroid translating along its minor axis 

As a third example, we determine the motion of a prolate spheroid translating along its minor axis. We consider 
the same prolate spheroid as in the previous example, and compute the swimming velocity U = Ux. when two jets 
(with fluxes ±q) are placed at x = ±6x. The resistance problem has a tractable solution and representation using the 
singularity methods described by Chwang and Wu [46] . In that work it was shown that a prolate spheroid translating 
along its minor axis with velocity U — t/x generates flow and pressure fields given by 



u — Ux — aBiii. — ax [ — — ) z — axrB^r + V ■^ fix 

i?2 ill 



-Ri- 



z + c 



R2 + BI 



p — 2/iQ;- 



X z — c z -\- c 



i?2 Ri 

where c = vt?— &^, f = \/x^ + U^, f = (^^x + yy)/r, i?i = y^z + cp + r^, R2 ~ \/(z — c)^ + r^, 



Bi = log 



i?2 



(z-c) 



Ri-{z + c) 



B, = 



I I z + c z — c 



Ri 



i?2 



and finally. 



2(3e^ 
1-e 



= Ue^ [e + (3e2-l)tanh"^(e)]" 



(23) 
(24) 

(25) 

(26) 



with e = c/a the eccentricity. From the flow and pressure fields the stress tensor d' may be computed without much 
difficulty. Integrating the stress over the body surface results in the expression 



F = 67r^Ct/x, 



(27) 



where 



C = ^e3 [g ^ (3g2 _ ^) tanh-i(e)] "' , 



(28) 



We need only determine the fluid stress in the resistance problem at the jet locations x = ±&x in order to compute 
the velocity in the swimming problem. After some algebra, we flnd that (with n = ±x at the jet locations). 



(n • 5- • n) (x = ±fe x) = — p + 2^ (n • E • fi) 

4 (2e{e'^-2) e(eC7 - a) + (a - Se^a) tanh"^(e) 

^ n 2 ^ ^^' 

avl — e. 



a^/^^^e2 [e + (Se^ - I) tanh"^(e)] 



(29) 
(30) 



Inserting the above into Eq. (10) and simplifying, the swimming velocity of the jetting prolate ellipsoid is found to be 

U^-7^- (31) 

F. Example 4: A viscous pinwheel 

A sufficiently asymmetric body can be made to rotate by appropriate placement of the jetting nozzles. The rotational 
velocity of a prolate spheroid driven about its minor axis, a "viscous pinwheel," can be determined using Eq. (10). 
Consider the same prolate spheroid as described in the previous example, but with its major axis aligned with the x 
axis. The body volume is V{e) = 7ra'^(l — e^), with e = \/a? — IP- ja the eccentricity. Pure rotation about the z axis 
may be driven with outward flowing jets at Xi = (x, — r(x), 0) and X2 = (—a;, r(x), 0), with r{x) = \/\ — e^\/ a^ — x^ . 
In order to maintain internal volume conservation, suction jets (inward flowing jets) are placed symmetrically at 
X3 = (x,r(a:),0) and X4 = (— x, — r(x), 0), and we write (?!)(x) = g((5(x, xi) + (5(x, X2) — (5(x,X3) — (5(x, X4)). The 
flow at all four nozzles act with equal strength to drive the rotation. The corresponding unit normal vectors are 
fio = (— r'(±a;)x =F y)/-\/l + ^''(a;)^ at the outward flowing jets, and n^ = (— r'(±a;)x ± y)/-\/l + r'(x)2 at the inward 
flowing jets. 

The torque required to rotate the prolate spheroid about the z axis with angular velocity il (the resistance problem) 
is 

L ^ /6.Mae)3(l-e^)ri .^ ^3^^ 



3[e-(l-e2)tanh-^(e)]' 



(see Refs. [47, 48]). During this rotation, the pointwise fluid stress on the rotating solid, a, has an analytical expression 
which may be deduced using the singularity method of Chwang and Wu [46] (though it is unwieldy and not included 
here). Assuming that the nozzles are small, only the normal component of the traction at the nozzle locations are 
relevant in setting the rotational velocity in the swimming/jetting problem, H = fi z, and so we have 



3[e-(l-e2)tanh-^(e) 
167r^(ae)3(l-e2)fi ^ 



d e-^i-e-jtann \e)\ ^^ , 



Figure 3a shows the rotational velocity of a spheroid with aspect ratio hja = 1/2 (e « 0.866) as a function of the jet 
placement distance from the particle center along the major axis. Here il is normalized by the relative flux per body 
volume. The result confirms intuition: there is no rotation when the jets are placed symmetrically along either major 
or minor axis, and hence there exists an optimal jet placement for inducing a body rotation. Figure 3b shows the 
maximal angular speed, il* , which may be obtained by appropriate jet placement for bodies of eccentricity e, along 
with the optimal jet placement location x* (inset). As the body becomes spherical the optimal jet placement limits 
to x* I a = l/v2, though the rotation response becomes increasingly small. As the body becomes more elongated 
the optimal jet placement moves closer to the poles, where the resulting forces can produce the largest torque on the 
body; along with the vanishing torque on a rotating slender body as e — > 1, the maximal rotational velocity increases 
without bound. 



G. Hydrodynamic efficiency 

We have shown that the velocities of a jetting body are dependent only upon the normal tractions in corresponding 
resistance problems, through Eq. (10). The issue of hydrodynamic efficiency is rather more involved. We define 
here an efficiency used in many other works in low Reynolds number locomotion, and consider for simplicity only 
translational motion in what follows (see Ref. [3]). The efhciency is defined as a ratio comparing the rate of mechanical 
work done onto the fluid in the resistance problem to that done in the swimming problem when the bodies are moving 



V{e) 




n 1 



V{e) 



n* 4 



0.4 _ 0.6 

x/a 




FIG. 3. (a) Rotational velocity response to two jet-pairs positioned symmetrically at ±x, for a body of aspect ratio b/a = 
^/l — €?■ = 1/2. The rotational velocity is normalized by the relative flux per body volume, (b) Normalized maximal rotational 
velocity as a function of eccentricity e (with V(e) = 7ra"^(l — e^)), and corresponding nozzle placement position (inset). 



at the same velocity U, 



E^ 



U F 



/ u-(-cr-n)dS' 
JdD 



(34) 



The difficulty in computing the efficiency, compared to simply computing the swimming velocity, is seen plainly in 
Eq. (34). Instead of integrating the fluid velocity in the swimming problem against the stress in the resistance problem, 
as in Eq. (8), here we must integrate against the stress in the swimming problem which depends intricately upon the 
precise form of the jetting flow profile. A theoretical bound on the above measure of £ ~ 75% was found by Stone 
and Samuel when the surface deformations act tangentially to the surface, and the body is spherical [9] . Since f is a 
mechanical and not a thermodynamic efficiency, motion with f > 1 is theoretically possible, in particular when the 
no-slip condition does not hold everywhere upon the body surface and there are sources or sinks of fluid or surface 
material. For example, Leshansky et al. [17] showed that a slender treadmilling spheroidal body can locomote with 
arbitrarily large efficiency by continuously introducing and removing surface material at the poles. 



Since the body is force and torque free in the swimming problem, the efficiency still requires only knowledge of the 
force component normal to the porous surfaces, and we have 



£ = 



U F 



E 



(35) 



am 



(x) (-n • o- • n) dS 



As expected, the work done by the swimmer is the work done by the jets against the normal stresses in the fluid. Note 
that Eq. (35) is only the external efficiency, and it does not include, for example, the internal work done to create 
the jetting flows in the first place. In the remainder of the paper we focus on axisymmetric spheroids, for which the 
swimming velocity and efficiency may be determined analytically, and we maximize the swimming efficiency through 
numerical optimization. 
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III. SPHERICAL BODY SHAPE 



General solution 



In the previous section we showed that the velocity of a spherical body is maximized by placing two small jets at 
the poles along the axis of locomotion. We now present an analysis of the dynamics and efficiency of a spherical body 
of radius a with an arbitrary axisymmetric fluid-jetting profile. We consider only jetting profiles which are fore/aft 
asymmetric and therefore volume conserving. We also assume the jets, and the resulting flow, to be axisymmetric. 
Swimming spheres with arbitrary velocity boundary conditions have been studied in a more general setting by Lighthill 
[49] and Blake [50]. Magar and Pedley [51], and Ishikawa et al. [52] have recently considered the behavior of swimming 
spheres with only tangential surface distortions, so-called squirmers, also using a similar analysis to that presented 
below. 

We set xi = z, X2 = — z, Sg = Sg = Ss, and q = —qi = (72 as before. Hence, the spherical body translates along 
the z direction and the Stokes equations may be solved in axisymmetric spherical coordinates (r, 0), with r G [a, 00) 
the radial distance and 9 G [0,7r] the polar angle. For notational convenience we again deflne ( = cos{9). The 
porous surfaces are taken to have spherical cap heights of lengths ea, as illustrated in Fig. 4, with e G [0, 1]. Hence, 
Ss = 2-Ec?e, and the polar angle illustrated in Fig. 4 has (^0 = cos(6'o) = 1 — e. 

In order to determine the efflciency E we will solve the Stokes equations, Eqs. (l)-(2), with the general boundary 
conditions, Eqs. (3)-(4). Exploiting the linearity of the Stokes equations, we decompose the problem into two parts, 
setting il) — ip' + ipg. The first problem (for ijj') corresponds to axisymmetric rigid body motion in the swimming 
direction, while the second problem (for tps) corresponds to the flow generated by a jetting sphere which is flxcd in 
space centered at the origin. The separation of the problem into a rigid body motion and a jetting flow component 
is the same separation used to determine the swimming velocity in the previous section. 

The flrst problem has a well-known solution. The stream-function for flow due to a solid sphere translating through 
a quiescent fluid with velocity C/z is 



^' 



U 



r' 1 - C^ 



&-<') 



(36) 



with associated pressure and normal strain-rate tensor component (see [44]) 



^'^2a^^^' f.E-f = £;;, = 0. 



(37) 




FIG. 4. A fluid-jetting sphere of radius a. The porous surfaces have equivalent spherical cap heights of lengths e a, with 
e £ [0, 1]. The surface area of each porous cap is Ss = 27ra^e, and i^o = cos(So) — 1 ~ e. 
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In order to solve the second problem for ips, we write the fluid velocity in a fixed spherical coordinate system as 
u = Uj.r + ugeg, and define an axisymmctric stream- function such that 

1 dips 1 dips .„„^ 

r"^ sm(6^) oO rsm(O) or 

The use of the stream-function ensures that the incompressibility condition (Eq. 2) is automatically satisfied, and 
Eq. (1) becomes 

D^ {D^^s) = 0, (39) 

where 

The boundary conditions on the body surface are the jetting profile in the direction normal to the surface, and zero 
tangential surface velocity, 

Ur{r = a,9) = c^iO, ugir = a,e)=0. (41) 

The general solution to Eq. (39), for flows which decay in the far-field and have bounded tangential velocities at the 
poles, may be written as 



^s-y: 






n = 2 



GniC), (42) 



with ^„(C) the n Gegenbauer function of the first kind [44, 53]. 



The Gegenbauer functions are related to the Legendre polynomials PniC) via 

5n{C) = ^^^^ (43) 

for n > 2. Under this representation the radial and tangential components of velocity are written as 

-r = ^f:(^ + ^)Pn-AC), (44) 



CO 






-^--E[(-^)^ + (--^)^)^=,- (45) 



We denote the inner product of the radial velocity boundary condition on the body surface and the n Legendre 
polynomial by c„: 

Cn = {Ur,Pn)=J cPiC) Pn{0 d^. (46) 

Thus, by successive inner products of Eqs. (44-45) we find 

(2n-l)(n-3) „+i (2n-l)(n-l) „_i 
An ^ a ^ c„_i, Bn = a c„_i. (47) 
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The corresponding pressure field, which may be found by integrating Eq. (1), is 



n=2 
1 



2(2n-3)B„ 

-^Pr^-liO, 



and the symmetric rate-of-strain tensor E = - (Vu + Vu"^) has an f f component 



Err = f • E • f = (?r 



n=2 ^ 



(48) 



(49) 



B. Efficiency and optimization of an entirely porous spliere 



The swimming efficiency £ from Eq. (34) may now be determined by combining the two problems, ip = ip' + ips- 
The rate of mechanical work performed on the fluid in the full swimming problem is 



$ 



/ U-{-(T-h)dS^ -2TTa^ I 0(C) [-{p + p') + 2flErr] d(, 
JdD J-1 



(50) 



where we have used that the swimming body imparts no force or torque onto the fluid. Inserting the expressions from 
above, we have 



-{p + p') + 2^,Err — M /, 



2{n + l)^+(2n + 2--\^ 



Pn-iiO-^fiUC, 



(51) 



and after further simplification and a shift in the summation we find 

3 " 



$ = 7ra^y^(2n + 1) 



(2n + l) + 



n+1 



4 + 37:^iaU C0(C)dC 



(52) 



The swimming velocity, from Eq. (10), and related towing force are 

^1 ^{^){i-H^))ds = -\j CHOdC, 



u = - 



F = 671 fi aU = — 37r/x a j C ((>{(,) dC,- 



(53) 

(54) 



The sphere moves with a non-zero velocity if and only if the jetting profile (j>{() contains the first Legendre polynomial 
mode. Moreover, this is the only component of the jetting velocity that contributes to the swimming speed. The 
hydrodynamic efficiency, on the other hand, does depend upon the full nature of the jetting profile 0(C) j and we find 



£ 



1 ^2 

C^iOdC 



sfJ^CcPiOdA +^f2{2n + l)(2n+l + -^\ ij /" 



{04>{C)dC 



n=2 



(55) 



As a simple example consider a sphere which is entirely porous (e = 1), which is propelled with the jetting profile 

0(C) = ^Pi(C)- Ac (56) 
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FIG. 5. The optimal jetting sphere (cross-section), with (j){() — Pi{i^)/Ss- (a) Strcamhnes for a jetting sphere which is fixed at 
the origin, u(r — a,9) = <j){{^)n. (b) Streamhnes for a swimming sphere, u(r = a,6) — Uz + (j>{Qn. (c) Distribution of surface 
fluid velocity for the fixed body, as in (a), (d) Distribution of surface velocity for the swimming body, as in (b). 



The jetting profile satisfies the flux constraints 



dD 



(x) dS = 2na^ [ </>(C) dC ^ 0, [ 0(x) dS ^^ q, [ 0(x) dS == -q. 



(57) 



Streamlines generated by the jetting body held fixed at the origin are shown in Fig. 5a (for q < 0), and in Fig. 5b 
when the body is swimming freely through the fiuid. The dipolar fluid structure in the latter is evident. The 
distribution of fiuid velocities at the body surface in each case are shown in Figs. 5c-d. The corresponding swimming 
velocity and hydrodynaniic efficiency are 



U 



Stto-^ 



£ 



1 



12.5% 



(58) 



The swimming speed is smaller than the largest possible swimming speed as discussed in the previous section by a 
factor of 2/3. However, it is simple to see that this jetting profile is in fact the most hydrodyamically efficient jetting 
profile possible. Every component of 0(C) which has a non-zero inner product with the n*'* Legendre polynomial, for 
n > 2, contributes positively to the denominator of £ in Eq. (55), without increasing the swimming velocity. Hence, 
the optimal surface jet pattern will only have the first Legendre polynomial component. Moreover, the decrease in 
the efficiency with the inclusion of any other Legendre modes in the jetting profile is dramatic, as can be seen by 
the form of Eq. (55). The higher Legendre modes do not contribute to the swimming velocity, but introduce large 
wave-number variations in the fluid flow, corresponding to an increased viscous dissipation of energy into the fluid. As 
an example, streamlines associated with the third Legendre mode are shown in Fig. 6, from which the corresponding 
increase in viscous dissipation may be intuited. 
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FIG. 6. Streamlines for flow generated by jetting profile 4>{0 oc PaiO- 



C. Efficiency and optimization of a partially porous sphere 

Organisms or manmade swimming devices utilizing non-inertial jet propulsion might be constrained to draw in and 
expel fluid from only some parts of the surface. Here we consider the efficiency and optimal jetting profiles of spherical 
bodies with e < 1 (see Fig. 4). 

The efficiency derived in the previous section, Eq. (55), is general and applies to such a partially porous body. 
We first note that the swimming efficiency is zero if the surface fluid velocity is discontinuous in (. Physically, this 
corresponds to solutions with non-integrable stress singularities at the edge of the pore (the Legendre modes needed to 
represent (j){Q do not decay at a sufficient rate to give a finite value of the efficiency). A discontinuous jetting profile 
is considered in Appendix A. The Legendre modes decay as c„ ^ n^'^", and as expected the sum in the denominator 
of Eq. (55) diverges. 

To see the general scaling of the efficiency with the porous surface area, consider the following form for the jetting 
profile which has (j){() = at the porous surface edges (so that the surface fluid velocity is continuous on (^ = cos(6') € 
f-1,11): 



'/'(C) 



Tra^ e 



1- 



(1-0 



X{CG[l-e,l]} 



- 1- 



(1 + 



X{Ce[-i,-i+e]} 



(59) 



where x is the indicator function, which would correspond to the quadratic profile of pressure driven flow through a 
single open pore as e — > (see Ref . [44] ) . The corresponding swimming velocity is 



U 



C^iC) dC 



-(3-e)g ^ 



(60) 



and we have 



'-^pnf (i_(i^)p„,c).c 



(61) 



The resulting efficiency f as a function of e is shown in Fig. 7a as a solid line. The efficiency is monotonically 
increasing with the size of the porous cap height, and limits to the optimal value for a spherical jetting body of 
£ — 1/8 = 12.5% as the body becomes entirely porous (e — >■ 1). Even for e = 0.2 (6*0 ~ 36°) the body achieves 
efficiencies on the same order as those mechanisms employed by biological organisms, £ ~ 3%. 

The optimal jetting profile cannot be determined as directly as in the entirely porous case, which was determined 
by simple inspection of Eq. (55). Instead, we seek numerically the optimal jetting profile by discretizing 0(C) at points 
corresponding to the Gaussian quadrature nodes of order M, and computing the efficiency via Eq. (55) keeping the 
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FIG. 7. (a) Efficiency as a function of e for the jetting profile <^((") — — 3- (1 — (1 — (^)/e) X{ce[i-e,i]} ^ (1 ^ (1 + 

C)/f)X{ce[-i,-i+e]} s^s a solid line, and the numerically determined optimal efficiency as a dashed line, (b) Optimal jet- 
ting profiles. Relative distance from the body indicates jetting fiow velocity at the surface. Individual profiles are scaled for 
presentation purposes, (c) Optimal jetting profiles, <j)*{0 ^s a function of (" for e = 0.5, e = 0.2, and e = 0.1. (d) Same as in (c), 
but scaled. The optimal flow profile through an increasingly small porous surface area appears to limit to the noted function 
and is shown as a dashed line. 



first N terms in the summation. The profile is subject to the constraints of Eq. (57), and (/)(C) = for |C| ^ 1 — £• 
The optimal profile is selected using an SQP quasi-Newton line search method built into the MATLAB optimization 
toolbox. The number of discritization nodes M and summation terms TV are increased until there is no discernible 
variation in the result. 

The optimal jetting profiles so obtained, denoted by 0*(C), are shown in Figs. 7b-d, and correspond to the efficiencies 
shown as a dashed line in Fig. 7a. The entirely porous case, e = 1, returns the expected result, a profile which scales 
linearly with C. Under the scaling indicated by the axes in Fig. 7d, it appears that the optimal jetting distribution 
through small pores at the poles (e — > 0) is limiting to the form 0(C) — 3(7/(47ra^e)Y/r^^^T^-(0A) which is shown as 
a dashed line. The maximal efficiency shown in Fig. 7a does not deviate dramatically from the example of previous 
consideration. 



IV. PROLATE BODY SHAPE 



The efficiency of swimming at low Reynolds number is known to depend significantly upon body shape. The drag 
on a solid body of a given volume is minimized when the shape is approximately a prolate spheroid with aspect ratio 
« 1/2, but with conical endpoints [54, 55]. Swimming organisms can enjoy a decreased fluid drag by selecting a 
similarly streamlined shape. However, in a departure from the more common shape-drag relationship, the propulsive 
mechanism of a fluid-jetting body is such that the efficiency decreases as the body becomes more prolate, as we shall 
show. We begin by solving for the efficiency exactly for one particular jetting profile, then move on to consider a 
more general framework for studying the fluid-jetting dynamics and efficiency. Finally, as in the case of a spherical 
body shape, we will determine numerically the optimal jetting profile for both entirely and partially porous jetting 
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surfaces. 



A. An exact solution 

We proceed just as in the spherical case, but the Stokes equations are now solved in prolate spheroidal coordinates 
[44]. The appropriate coordinate system is given by the conformal mapping 

z + ir ^ ccosh{£_ + i 9), (62) 

and we set 

T = cosh(0, C = cos(6') (63) 

for clarity. The surface defined by r = tq > 1 is a confocal spheroid with foci ±cz, and we have r G [1, oo), ^ G [—1,1]. 
For a spheroid with major and minor axis lengths 2a and 26, respectively, we have tq = a/c = a/\/a'^ — b^, the inverse 
eccentricity as in §11. The Stokes stream-function satisfies the biharmonic equation 

D'^ij = 0, (64) 

where 



D^^ 



1 



C2 (t2 - C2) 

Writing the velocity in the new coordinates, u = UrT + u^^, where 



(r2 - 1) drr + (1 - C') ^cc . (65) 





r 




- 1 


CV^ 


- 1 




-c% 



V^^^e Vr' - c 



(66) 

(67) 



we set 



1 dip 1 dip 

Ut = :— -, Ur = :^:— , (68) 

CV(^' - C") (^' - 1) ^^ cV(t2-C2)(1-C')9t' 

so that the incompressibility condition is automatically satisfied. Also of use are the relations d^ = (t^ — l) dr 

and dr, = - {l - Cy^^ di^. The surface area element is dS = J{C)dCde, with J(C) = cV(to - l)(^o - C^) = 
by/a^-(a^-b^)C. 

An exact solution can be derived for the particular jetting profile 

which limits to the linear profile as studied in the spherical case, Eq. (56), as 6 — t- a. Once again we split the problem 
into two parts, ip = ip' + tps- The first problem for ip' is again that of rigid body motion along the axis of symmetry 
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with swimming velocity Uz. Here also the solution is known [44], and we have 



^' = ~^(1-C^)(r^-1) 



(-0^ + 1) / {rS - 


- 1) log 


V + l' 
r-1 


2r 


t2-1 


(^o' + 1)/(to'- 


-l)log 


ro + 1 

ro-lj 


2to 


t2 - 1 



(70) 



The corresponding pressure field may be found by integrating Eq. (1) in the appropriate coordinates, leading to 



M 



dp 

dr c(t^ — 
dp M 



dC c(l - e) 






yielding, upon integration, 



p'(t,C)- 



2^iUC 



KT2-C2)[(T2 + l)cOth-l(ro)-To] 

The TT component of the symmetric rate-of-strain tensor is 



(71) 
(72) 

(73) 



S' = -f • E' • f = 



7^2^71; Qu^ C^/l^ uc_ 



c(t2-C2)1/2 dr c(t2-C2)^/^ 



(74) 



Inserting Eq. (68) and simplifying, we find on the surface t = tq that _E^^(to, C) = 0. The pressure above corresponds 
to the stress used in the example of a prolate spheroid translating along its axis of symmetry in §11. As shown in 
that example, the stress on the body during rigid body motion is sufficient information for determining the swimming 
speed in the full swimming problem. 



From Eq. (10), we have as in the example of §11 (replacing U with U in p' and -E^^), 



^2 /•! 



^^-e^L^W^-^' + ^'^^-^^^-Sy-r 



0(C) i-p) 



{ri-l){Ti-e)dC, 



(75) 



leading to the general expression 



U ^-'- 



GiCmOdC, G(C) = C^ 






(76) 



As 6 — > a, we recover G(C) — C bs found for the spherical body. Inserting the jetting profile from Eq. (69), we finally 
obtain 



U 



gr^[rocoth ^(tq) - 1] 



(77) 



The velocity limits to that of the spherical body as tq ^ oo, and increases without bound (for fixed q) as the body 
becomes infinitely slender. 

In order to determine the efficiency, we must now solve the second problem (for tps), that of a jetting spheroid fixed 
in space at the origin. 



D*i>s = 0, 



(78) 
(79) 
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The form of the velocities on the body surface suggest the fohowing ansatz for the stream function, 

^^-(1-C')5(t), (80) 

and after some manipulations we arrive at the solution, 

W, ^ g (I .2N [^+(^'-l)'^Qth-^(r)-2rroCoth-^(ro)] 

27r^ ^ ^ (r2 + l)coth-i(To)-ro ^ ' 

The corresponding surface pressure is, upon integration of Eqs. (71)-(72), 

/ .X 2gCA^[7'oCOth~^(To)- 1] 

c^i- (tq - C^) [(^0 + 1) coth (ro) - ToJ 
and the normal component of the traction is (from Eq. (74)), 

£„ . ^i^"^' + i[-^) ^ ,83) 

We now have the required information for determining the hydrodynamic efficiency in the full swimming problem. 
The work done on the fluid by the swimming body is 

$ = -271 / 0(C) (-P + 2iiE,, - p')J{C) dC (84) 

(85) 

Meanwhile, using F = 67r/ii?U with R as defined in Eq. (19), we find the efficiency 

^^2K-l)(.„,otl.-(r„)-l)' ^^^ 

[to - (tq^ + 1) coth \to)] 

The efficiency is found to decrease monotonically to zero as the body becomes more slender, as shown in Fig. 9a as 
a dashed line. In the spherical limit (as 6 ^ a) we recover £ -^ 12.5% as expected. 

As previously noted, it is unusual that a moving body's hydrodynamic efficiency decreases monotonically with its 
slenderness. The reader may have already guessed this result, however, given the nature of the propulsive mechanism. 
The jetting flow through the body is constrained to act locally in the direction normal to the body surface. As the 
body becomes more prolate in shape, the nozzles on the surface become oriented in a direction more perpendicular 
to the direction of motion, and become much less useful for propulsion. From this observation, it is thus expected 
that the efficiency will increase as the body becomes more oblate in shape, even though the surface area presented to 
the fluid is increasing in that case. This will be the subject under consideration in §V. First, though, we develop a 
general framework for studying prolate bodies with arbitrary jetting profiles, and determine numerically the optimal 
actuation for both entirely and partially porous body surfaces. 



B. Efficiency and optimization of an entirely porous body 

The dynamics and efficiency of a prolate jetting body for an arbitrary jetting profile, and for partially porous 
surfaces, requires a different approach to solving the jetting problem for ips- The general solution to Eqs. (78)-(79) 
may be expressed using spheroidal harmonic functions i?„ and Sn, themselves expressed in terms of the Legendre 



271 y ^{C)i-p+^^^Err-p')J{C)dC 






c% li (r2 - ef (r2 - 1) - c3^ 


TO- 


-(T2 + l)cOth-HTo)" 
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FIG. 8. (a) The optimal geometry-normalized jetting profiles for a sphere (dashed line) and a prolate spheroidal body (solid 
line), (b) The same, for a sphere (dashed line) and an oblate spheroidal body (solid line). 



polynomials of the first and second kind, P„ and Qn respectively (see [56]): 



4'5 = Ai i?i(r)i?i(C) + ^ [Sn + S„ (t2 + C')] i?„(C)5„(r) 



(87) 



where 



R^{x) = (1 - x^)P'^{x) = -nxP„(a;) +nP„_i(a;), 

Snix) = (1 - x'^)Q'„{x) = -nxQnix) + nQn-i{x), 

dX 



Qn{x) = Pn{x) 



P„(A)2(A2 - 1) 



(89) 
(90) 



As discussed in more detail in Appendix B, the boundary conditions are used to determine the coefficients i3„ and 
Bn- Subsequently, the pressure is found through integration of Eqs. (71)-(72), and for arbitrary ((){(,) we find that 



ErriroX) 



ro(l + C^-2rg) 
c(r2-e)3/2(r2- 1)1/2 



0(C). 



(91) 



Meanwhile, the solution to the rigid body problem for ip' (Eq. 70) does not depend on the jetting profile outside 
of the dependence of U on 0(C), and the results of the previous section carry over here. For a given jetting profile 
(/)((■), the first N coefficients B^ and B^ are determined numerically and the efficiency is then found using a Gaussian 
quadrature (see Appendix B). 

Given the framework described above, the optimal jetting profile was selected by a similar SQP quasi-Newton line 
search method as the one used for optimizing the spherical swimmer. The optimal profile '/>*(C) so found (normalized 
by the geometric scaling J{C,)/q) is shown for a body of aspect ratio h/a — 0.3 in Fig. 8a, along with the optimal 
profile for the spherical case {h/a =1). As the body becomes more slender in shape, the only jets which contribute 
appreciably to the swimming velocity, via Eq. (76), are those near the fore and aft poles. However, along with 
increased jetting velocities near the poles comes dramatically increasing viscous stresses there, and hence more work 
is done on the fiuid. Since J(C) ^ eyl — C/^ for e = h/a <C 1, the optimal profile shown in Fig. 8a requires a much 
more significant placement of jetting nozzles near the fore and aft poles, with a slight reduction of flow strength at 

C = ±i. 

The optimal efficiencies for a range of aspect ratios are shown in Fig. 9a as filled circles, along with the exact 
efficiency calculated for the particular jetting profile 0(C) — qC, / [tt J {(/)] from Eq. (86) as a dashed line. (The other 
plotted values correspond to partially porous bodies, as discussed in the following section). While the optimal profile 
differs from that as studied analytically in the previous section, the efficiency increases only very slightly for each 
aspect ratio, by approximately 10""^. The integrated form of the efficiency measure appears to be insensitive to 
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FIG. 9. (a) Optimal swimming efficiency as a function of aspect ratio for a selection of nozzle sizes, limiting to £* = 1/8 
as the body becomes spherical and entirely porous. The efficiency given by the jetting profile from the previous section 
(0(C) = qC / {''^ J {()) ^ e = 1) is shown as a dashed line, (b) Maximum efficiency for a given aspect ratio b/a and nozzle size e. 
Contours are shown at efficiency value multiples of 0.01. 



the precise form of the jetting flow profile for prolate jetting bodies, provided the general behavior of approximate 
geometry-normalized linear growth in (^. 

C. Efficiency and optimization of a partially porous body 

Just as in the spherical case, we can ask about the optimal jetting profile and maximal efficiency when the swimming 
body's surface is only partially porous. Here we assume that the body is porous near the fore and aft poles in the 
regions |C| > 1 — e (or \z\ > (1 — e)a), and satisfies a no-slip condition along |CI < 1 — £■ The parameter e is the 
analogue of the spherical cap height parameter from Fig. 4, and the vertical measure of the porous cap height is again 
ea. The optimization is performed numerically as before, but including the constraint that 0(C) — for |C| < 1 — £■ 

The maximum efficiency, £* , is shown for a range of aspect ratios b/a and dimensionless cap-lengths e in Figs. 9a-b. 
For a given e the efficiency decreases as the body becomes more slender for the reasons already discussed. For a given 
aspect ratio, the maximal efficiency must not increase as the body becomes less porous, since the optimal profile for 
a more porous body could have 0(C) ~ wherever necessary. The efficiency becomes less sensitive to the porous 
surface size as the body becomes more prolate in shape, even though the jetting profile changes rather significantly as 
e is decreased. Compared to the efficiencies of many other propulsive mechanisms at low Reynolds number, including 
those exploited by nature, the fluid jetting body compares favorably for all but the most slender of bodies, or the 
smallest of porous cap lengths. 

The optimal jetting profiles are illustrated for a body with b/a = 0.2 for e = 1, e = 0.5, and e = 0.2 in Fig. 10a. 
For each value of e, the distance of the curve from the body surface corresponds to the jetting flow speed there. (The 
distance is scaled differently in each case for clarity of presentation.) The profile varies slowly in C in the fully porous 
case, and much more dramatically in the e — 0.2 case. However, all three profiles correspond very nearly to f « 1% 
as indicated in Fig. 9a. 

V. OBLATE BODY SHAPE 



Generally, efficient swimming in viscous fiuids requires a streamlined body shape, and drag anisotropy is frequently 
utilized for propulsion. We have seen in the previous section that the hydrodynamic efficiency of a jetting body 
decreases as the body shape becomes more prolate, due to the reduction in the presented surface area in the direction 
of motion. The jet nozzles becomes less aligned with the swimming direction, and only the nozzles near the poles 
may be used for propulsion, adding significantly to the hydrodynamic work done on the fluid in the process. We 
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FIG. 10. The optimal jetting profiles are illustrated for a (a) prolate spheroid with hja — 0.2, (b) sphere, and (c) oblate 
spheroid with a/b = 8. For each value of e, the distance of the curve from the body surface corresponds to the jetting flow 
speed there. (The distance is scaled differently in each case for clarity of presentation.) 



now explore bodies in the shape of oblate spheroids, and determine the consequences on swimming efficiency. We 
proceed just as in §IV, benefitting from a simple correspondence between the prolate spheroidal and oblate spheroidal 
coordinate systems. 



An exact solution 



The results from the previous section may be used to determine the swimming velocity and efficiency for an oblate 
spheroid which is moving along the direction of its symmetric axis. The prolate spheroidal coordinate system is 
transformed to an oblate spheroidal coordinate system by replacing c by ic, r by — iA, and reversing the roles of a and 
b so that a > fe in both cases. The surface area element becomes dS — J(C) d(^ d9, with J(C) = c^ ■^/(Aq + 1)(Aq + C^) — 

The exact solution in the prolate case is simply converted to the oblate spherical coordinate system, and we assume 
the analogous jetting profile (f>{C) = qC/[TrJ{C)]- The swimming velocity is then 



U = 



-q[l-XoCot-^{\a)] 



(92) 



As the body becomes spherical (Aq — > oo) the swimming speed limits to the expected value U = —q/{i'Ka?). If instead 
we consider the limiting case of a flat-plate moving along its axis of symmetry (6 — > a, Ao — > 0), we find that the 
swimming speed becomes U — —q/{iTa^). This result may not be surprising, given that as Aq — >■ the jetting profile 
becomes (/>{(,) = q/{TTa^)H{(), where -ff(C) is the Heaviside function with a jump at C = 0. The flow speed is constant 
everywhere on the body surface outside the disk edge C = 0, and the swimming speed is precisely that opposing this 
jetting flow speed. 

The efficiency of an oblate spheroidal jetting body with this jetting profile becomes, transforming Eq. (86), 



£ = 



2(Ag + l)[AoCot-i(Ao)-l]^ 
[(A2-l)cot-i(Ao)-Ao]' 



(93) 



In the spherical limit the efficiency becomes £ — 12.5% as expected, and as the body becomes a flat plate we find 



lim £ — 



81%. 



(94) 
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FIG. 11. (a) Optimal swimming efficiency as a function of aspect ratio for a selection of nozzle sizes, limiting to £* = 1/8 
as the body becomes spherical and entirely porous. The efficiency given by the jetting profile from the previous section 
{4>{0 = qC / {''^ J iO) ^ e = 1) is shown as a dashed line, (b) Maximum efficiency for a given aspect ratio b/a and nozzle size e. 
Contours are shown at efficiency value multiples of 0.1. 



Even without further optimization of the jetting velocity profile, the jetting body outperforms such ubiquitous propul- 
sive mechanisms as flagellar undulations by an order of magnitude in this measure. 

B. Efficiency and optimization of an oblate jetting body 



The jetting profile is now optimized numerically as in the previous section. Figure 11a shows the maximal efficiency 
found as a function of the aspect ratio a/b for an entirely porous body (e = 1), along with the efficiency from the 
exact values from the previous section for the particular jetting profile considered there (as a dashed line). Unlike 
for prolate bodies, here there is a dramatic increase in the efficiency for very oblate bodies upon optimization. The 
jetting profile used to find an exact solution in the previous section limits to a Heaviside function as the body becomes 
a flat-plate. A consequence of this dramatic jump in the surface jetting velocity is a high wavenumber dissipation, 
with its associated deleterious effect on the swimming efficiency. Perhaps, then, the optimal jetting profile should 
have a tapered magnitude near the curve ^ = 0. This suspicion is confirmed numerically, as shown in Fig. 8b for a 
jetting body with aspect ratio a/b = 8. Outside a region near C = the optimal profile (normalized by J{C/)/q) is 
approximately linear in C,. 

Numerically, we find that the computational results for increasingly oblate bodies falls increasingly close to the 
approximate jetting profile 
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(95) 



We may insert Eq. (95) into the framework developed above and determine the corresponding swimming efficiency for 
any range of aspect ratios. The result is shown in Fig. 12. The efficiency corresponding to the above jetting profile 
is significantly larger than that of the previously considered analytical case, appearing to slowly limit as the body 
becomes a flat-plate to a value of f sa 162%, which is well above one. This value is also confirmed by our optimization 
results. For oblate-shaped bodies, the swimming efficiency is therefore quite sensitive to the (geometry-normalized) 
form of the jetting profile. 

Finally, the partially porous case is considered for oblate spheroidal bodies. The maximal efficiencies are shown in 
Figs, lla-b for a range of aspect ratios a/b and cap distances e (where the surface is porous again on |C| > 1 — e, or 
\z\ > (1 — e)b). The optimal profiles are illustrated in Fig. 10c for e = 1 (entirely porous) and e = 0.5. As in the 
previous considerations the efficiency can only decrease with decreasing e. However, for a given porous surface area 
we find that the efficiency does not vary dramatically as the body becomes more oblate. When the entire surface is 
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FIG. 12. The efficiency corresponding to two assumed jetting profiles for a large range of aspect ratios. Dashed line: analytical 
profile studied in §VA. Solid line: approximate jetting profile leading to the majcimum swimming efficiency (Eq. 95), which 
yields an efficiency tending towards approximately £ ~ 162% as the body becomes a flat plate. 
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porous we have already seen that tapering the jet strength near C — can have a sizable impact on the swimming 
efficiency. If the body is only partially porous, the jetting profile near the edges is already zero, and the primary 
consequence of increasing the oblate aspect ratio must be a variation in the work done in the resistance problem. Since 
the drag on an oblate spheroidal body changes very little with the body aspect ratio, the slowly varying efficiency 
profiles might not be particularly surprising. 



VI. DISCUSSION 

A. Summary of our results 

In this paper we have shown that a body can swim with remarkable efficiency by drawing in and expelling fluid 
at zero Reynolds number using non-inertial jet propulsion. Optimization of the jetting velocity profile at the surface 
yielded (external) swimming efficiencies oi £ — i2.5% for a fully-porous spherical body, and as much as £ « 162% as 
the body becomes very oblate in shape. In comparison to other common mechanisms utilized at low Reynolds numbers 
such as flagellar undulations, non-inertial jet propulsion thus presents an improvement of two orders of magnitude 
[6, 7]. 



B. A locomotion mechanism for all Reynolds numbers 



An interesting aspect of jet propulsion, unlike such mechanisms as an undulating flagellum, is that a fluid-jetting 
body can self-propel at all Reynolds numbers, and it does so by exploiting distinct physical regimes. For large Reynolds 
numbers the motion of bodies generated by the expulsion of fluid jets has generated a vast literature, both in the 
exploration of jet propulsion in nature and for engineering purposes (see for example Refs. [36, 38-42]). Unlike at 
zero Reynolds number, where the swimming efficiency is independent of the fluid flux q, greater thrust and efficiency 
may be achieved at higher Reynolds numbers by tuning properly the vortex ring formation in a pulsatile jetting 
locomotion [42]. At low Reynolds numbers the non-inertial jet propulsion utilizes the generation of viscous stresses 
in the surrounding fluid for propulsion, while at higher Reynolds numbers the mechanism shifts to the transmission 
of momentum into the fluid opposite the body motion. In contrast, the locomotion of a rotating helix, for example, 
while effective in a Stokesian realm, is all but useless at higher Reynolds numbers. 
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C. The locomotion of Synechococcus 

A strain of cyanobacteria, Synechococcus, swims in fluid, but does so without the use of flageUa [9, 33-35]. In 
fact, their motihty mechanism is stiU an open problem in biophysics. Stone and Samuel [9] proposed a model where 
Synechococcus was assumed to locomote by compressive surface distortions, i.e. the passage of traveling waves, tan- 
gentially along the body surface. Here we have shown that another mechanism of motility, one which includes a 
surface flow which acts in the direction normal to the body, can also provide an efficient locomotion. The motility 
of slime-extruding organisms such as cyanobacteria and myxobacteria (see Fig. 1) appears experimentally to depend 
significantly on the non-Newtonian rheology of the extruded slime, and the presence of a substrate. However, Syne- 
chococcus is known to swim absent the presence of a substrate, without changing shape, and without any observable 
external organelles. 

Given the recent observations of porous fluid extrusion in the biological community [27-32] , the work presented here 
might suggest an alternative motility mechanism for Synechococcus, that of non-inertial jet propulsion. To this end, 
let us approximate the swimming speed and efficiency of Synechococcus under the assumption that it utilizes fluid 
extrusion as a propulsive mechanism. The organism shape is approximately that of a prolate spheroid with aspect 
ratio b/a sa 1/2, with cell length 2a w 10^'* cm. Wolgemuth et al. [31] note that a related organism has nozzles 
everywhere on the surface, but Fig. 9 indicates that the efficiency is not exceedingly dependent on the porous surface 
size. For the wide range 0.2 < e < 1 we recover an optimal swimming efficiency of 2% — 5%, well within the range 
of efficiencies reached by other microorganisms, if not slightly larger. To compute the swimming speed we need the 
flux q. To approximate the value of q we use the estimated figures of Wolgemuth et al. [31] for a related organism 
M. xanthus. Namely, that there are approximately N = 500 nozzles around the circumference of each end of the cell, 
each nozzle has cross-sectional area A w 3 x 10^^"^ cm^, and the fluid exit velocity \s u k, 10^'"' cm/s. This yields a 
fluid flux at each end of the cell of g = NuA « 1.5 x 10~^^ cm'^/s. From Eq. (77), and assuming that jet propulsion is 
indeed the main locomotion mechanism, the swimming velocity would be approximately \U\ ~ Q.2q/a? = 1.2 x lO"'' 
cm/s = .0012 /im/s. This swimming speed is many orders of magnitude smaller than the observed swimming speed 
of Synechococcus, which is closer to 10 iira/s. In order for the swimming velocity to match that of Synechococcus, the 
flux would have to be q « 10^^^ cm^/s, which would require a volume equivalent to the internal body volume to be 
expelled in an exceedingly fast time, t « 10~^s. We can therefore conclude that the expulsion of a Newtonian fluid 
is unlikely to be the propulsive mechanism utilized by Synechococcus, giving support to other theories such as the 
passage of traveling surface waves of tangential displacements and of small amplitude [9] . 

D. Directions for future work 

Exciting directions which this line of inquiry may take include optimization and control for jetting motions in three 
dimensions. For example, for a body with a known generalized resistance matrix, how best to utilize a finite collection 
of surface jets to move from one point to another in space? Separately, interaction dynamics may yield unexpected 
dynamics for low Reynolds number swimming behavior given the possibility of drawing in fluid just expelled by other 
swimmers. Finally, in the physical realization of such a swimmer, the internal mechanisms for driving the flow both 
into and out of the body must be addressed. A more detailed efficiency analysis would then be of great interest, where 
internal costs are considered in addition to the hydrodynamic work performed on the external fluid. 

We have seen that a fluid-jetting body can swim with remarkable efflciency. The constraint that fluid flow in and 
out along the direction everywhere normal to the body surface led to a bound on the efficiency of 12.5% for a spherical 
body, and approximately 162% as the body becomes a flat plate. Imagine this constraint to be removed. In this 
case the body could draw in and expel fluid everywhere upon its surface precisely in the direction of motion, and 
the surrounding fluid would go undisturbed (see Fig. 13). Hence the work done externally on the surrounding fluid 
would be zero, and the motion would correspond to an infinite hydrodynamic efficiency. Generally, the muted fluid 
disturbance even with the normal-jetting constraint hints at another very intriguing aspect of this form of locomotion: 
the body does not strongly signal its presence at any distance as it swims through the fluid. The non-inertial jetting 
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FIG. 13. An arbitrarily shaped body does not disturb the fluid and swims with infinite hydrodynamic efficiency if the surface 
distortion velocity is constant and unidirectional. 



body as portrayed in Fig. 13 would be the perfect stealth swimmer. 
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Appendix A: Discontinuous jetting profile 

A discontinuous jetting profile results in zero swimming efficiency. To see this, consider a spherical body of radius 
a with the jetting profile 



<^(0 - ^:^ (X{ce[o,i]} - X{Ce[-i.o]}) 



A-KO?' 



where x is the indicator function. The swimming velocity is, from Eq. (11), 
Inner products against the Legendre polynomials give 
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Thus, the sum in the denominator of Eq. (55) diverges. Taking the first N Legendre modes to approximate 4'{C), 
the efficiency is made arbitrarily small, with £ — > in the limit of the discontinuous jetting profile (/'(C)- A stress 
singularity associated with the jump in surface velocity results in infinite work being done on the fluid. 



Appendix B: Framework for general solution for spheroidal bodies 



We present here the framework used to determine the behavior and efficiency of prolate spheroidal jetting bodies. 
The methodology for oblate spheroidal bodies is identical, under the mapping {t,c) — > (— zA,zc), and switching the 
roles of a and b (see Fig. 8). 

The general solution to Eqs. (78)-(79) may be expressed using spheroidal harmonic functions i?„ and Sm themselves 
expressed in terms of the Legendre polynomials of the first and second kind, P„ and Qn respectively, 



^S - Ai i?l(T)i?i(C) + J2 K + ^- (^' + C')] RniOSnir), , 



(Bl) 



(see §IV). From the Legendre differential equation, we have R'^ix) — —n{n + l)Pn{x) and S'j^{x) = —n{n + l)Qn{x). 
For a decaying fluid velocity as r — )■ oo we must have Ai ~ 0. If the jetting profile 0(C) is assumed to be odd about 
C = 0, then we may set i?„ = i3„ = for all n even. 

Due to the exponential decay of the terms Qn{To) and 5'„(to) as n — )• oo, it is helpful to define the normalized 
coefficients 



noting the tractable behavior of the ratio 
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Differentiating the expression Eq. (BI), the boundary conditions (Eqs. (79)) therefore yield 
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We determine the coefficients C„ and C„ numerically by enforcing the conditions above at a finite number of nodes 
in C € [—1, 1], corresponding to the Gaussian quadrature nodes of order M, and keeping the first N terms in each 
summation. The resulting linear system is inverted using simple Gaussian elimination. 

With the coefficients C„ and C„ in hand (and hence _B„ and -B„), we may determine the pressure and the tt 
component of the rate-of-strain tensor E. Associated with the stream function in Eq. (BI), Eq. (72) gives 
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which is integrated to within machine precision accuracy using Gaussian quadrature. The integration constant may 
be set to zero since the pressure must be an odd function about C = 0- 
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FIG. 14. Exact efficiency as a function of aspect ratio (with a denoting the longer semi-axis length in both the prolate and 
oblate cases) is shown as a dashed line. The computed efficiencies using M — 100 and A'' = 40 are shown as filled circles, 
showing excellent agreement over a wide range of aspect ratios. 



From Eq. (74), the desired rate-of-strain tensor component may be written as 

Ett{tq,0 
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This component of the stress is dependent upon the jetting profile in a straightforward manner. Namely, using the 
(arbitrary) jetting boundary conditions for ipg, we find that 
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Finally, the efhciency is determined by numerically integrating the rate of work done on the fluid. 
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and using the definition of Eq. (34), 
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The terms p', R, and U are just those as shown in Eqs. (73), (19), and (76), respectively. 

The numerical approach is compared to the exact solutions of Eqs. (86) and (93) in Fig. 14, using the jetting profile 
(/>{(,) = (^q/[TrJ{(^)]. Here we used AI = 100 and N — 40. The computed efficiencies show excellent agreement with 
the analytical results over a wide range of body aspect ratios. 
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